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VOLUME  INTERSECTION  OF  TWO  IDENTICAL 
RIGHT  CIRCULAR  CONES 

Ronald  G.  Hughes 
J.  Dean  Clamons 

INTRODUCTION 

The  successful  operation  of  an  instrument  system  depends  on  signal- 
to-noise  levels  as  well  as  sensitivity.  During  the  design  phase  of 
engineering  hardware,  a  calculation  is  required  to  predict  the  appli¬ 
cability,  within  a  particular  set  of  physical  conditions,  of  the  total 
instrument  package  to  the  phenomena  to  be  measured.  This  provides  the 
constraints  or  limits,  related  to  sensitivity  and  signal- to-noise  level, 
which  must  be  met  by  each  of  the  instrument's  system  modules,  and  these 
serve  as  the  goals  for  module  parameters  necessary  to  produce  a  func¬ 
tional  instrument. 

As  applied  to  the  doppler  current  meter,  a  calculation  to  deter¬ 
mine  the  volume  reverberation  signal  strength  is  made,  based  on  the 
assumption  used  with  sonar  frequencies,  that  the  reverberation  signal 
strength  is  directly  proportional  to  the  common  intersection  volume  of 
the  beam  patterns  of  a  transmitter-receiver  pair1  (see  Fig.  1).  Con¬ 
centric  beams  result  in  an  infinite  volume,  excluding  consideration  of 
acoustic  signal  attenuation,  but  introduce  the  problem  of  testing  the 
transducer  system  in  a  finite  container  since  the  reflection  from  the 
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walls  results  in  a  noise  signal  which  far  exceeds  the  volume  reverber¬ 
ation.  But  when  the  beams  of  the  transducer  pair  are  crossed,  the 
desired  signal  is  observed  in  the  output  of  the  signal  processing 
electronics. 

Following  the  above,  the  magnitude  of  the  signal  expected  in  the 
receiving  transducer  of  the  doppler  meter  configuration  is  obtained  as 
a  volume  reverberation.  This  signal  strength  is  assumed  directly  pro¬ 
portional  to  the  common  volume  which  is  viewed  by  the  receiving  trans¬ 
ducer  and  insonified  by  the  transmitter,  presupposing  suppression  of 
side  lobes.  The  predicted  value  of  the  backscattered  signal  is  calcu¬ 
lated  as  a  product  of  factors  related  to  this  common  volume,  the  volume 
scattering  function  model,  the  beam  pattern  spreading,  and  the  attenu¬ 
ation  loss. 

To  implement  this  approach  in  predicting  signal  levels,  the  volume 
intersection  of  the  doppler  meter  transducer  pair  is  calculated  with  an 
assumption  that  both  the  radiated  and  received  beam  patterns  are  conical 
and  identical. 

MATHEMATICAL  FORMULATION 

The  calculation  of  the  volume  intersection  of  two  right  circular 
cones  is  very  simple  in  principle,  but  involves  the  tedious  evaluation 
of  a  triple  integral.  Utilization  of  the  symmetry  of  the  problem  in 
cylindrical  coordinates  fails  since  the  area  formed  by  the  intersection 
of  a  cone  with  a  plane  not  perpendicular  to  the  symmetry  axis  is  ellip¬ 
tical  rather  than  circular. 


Further  restrictions  imposed  on  the  problem  follow: 

1.  The  vertices  of  the  two  cones  are  positioned  on  the  x-y  plane 
such  that  Z  ■  y  ■>  0. 

2.  Tie  vertices  are  separated  by  a  distance  D  >  0. 

Variables  include: 

9,  the  angle  of  intersection  of  the  symmetry  axes  of  the  two  beams; 
Ip  divergence  angle  of  transducer  beam  (apex  angle  of  the  cone); 

Al,  radius  of  the  transducers;  and 
S,  separation  of  the  centers  of  the  transducers. 

The  analytic  expression  for  a  cone  having  its  axis  of  symmetry 
congruent  with  the  Z  axis  of  an  orthogonal  Cartesian  coordinate  system 
is  given  by 

x3/a3  +  y3/b3  -  —  =  0  (1) 

ca 

For  the  right  circular  case  eq.  (1)  changes  from  (see  Fig.  2)  to 

x3  +  y3  -CaZa  =  0  (2) 

The  constant  C  appearing  in  equation  (2)  is  defined  in  terms  of  the 
divergence  angle,  ip,  specified  for  the  transducers,  as 

C  =  tan  Ip  / 2. 

The  expression  for  the  rotated  and  translated  cone  in  terms  of 
the  Cartesian  coordinate  system  in  Fig.  1  is  determined  by  means  of 
a  rotation  and  translation  operation  applied  to  equation  (2).  The 
rotation  operators  about  the  y  axis  for  angles  -  0i  and  +  0i  for  cones 
1  and  2  respectively  are: 


(cos  0  0  -sin  0  \ 

sir.  8  0  cos  0  ' 


Ry  ■  /  cos  8  0  sin  0  ’ 

3/1  1 

I  0  1  0 

\-sin  8  0  cos  0  . 

i  l. 


The  angle  0  in  the  expressions  for  the  rotation  operators  is  defined 
in  terms  of  angle  0  as  noted  in  Fig.  1:  0^  K  0  -  ff/2  for  cone  1,  and 

0  =  tr/2  -  0  for  cone  2.  These  substitutions  in  Ry  and  Ry  result  in: 
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Operating  on  equation  (2)  with  Ry^  and  translating  along  the  x  axis  by 
-D/2  yields  for  cone  1: 

((x  +  D/2)  sin0  -  Zcos0)a  +  yS  ■  tan3  iJ)/2  ((x  +  D/2)  (3) 

cos0  +  Zsin0) 


A  similar  application  of  Ry^  and  translation  of  D/2  on  the  axis  produces 


the  equation  for  cone  2  as 


Terror 
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( (x  -  D/2)  sin0  +  Zcos 0)2  +  y3  =  tan3  0/ 2  ((x  -  D/2)  (4) 

cos 0  -  Zsln0)2 


Because  of  the  symmetry  of  the  intersection  volume  with  respect  to  the 

x  =  0  and  y  =  0  planes,  it  is  necessary  to  consider  only  that  part  of 

volume  which  lies  in  the  first  (x  -  y)  quadrant.  The  volume  under 

consideration  can  be  further  segmented  into  Volume  V  and  y  as  shown 

a  b 

in  Fig.  1.  Simplification  of  integration  results  in  particular  from 

setting  the  limits  as  V  is  bounded  by  cone  1  and  the  planes  x  =  0, 

y  *  0,  Z  =  Z  and  is  bounded  by  cone  2  and  the  planes  x  =  0,  y  *  0, 

Z  =  Z  .  The  integral  in  Cartesian  coordinates  for  Va  is  found  to  be: 
a 


V 

a 

4 


Z/tan(0  -  0/2)  -D/2 
dx 


X 


(tan2  0/2  ((D/2  +  x)  cos0  +Zsin0)3 

Zcos0)2  ) 


-  ((x  +  D/2)  Sin0  - 
i/a 


Integrating  this  expression  with  respect  to  the  y  variable  and  redefin 
ing  the  terms  in  the  resulting  integral  as 
A  =  sins0  tan2  0/2  -  cos20 
B  =  2  sin0  cos0  (tan2  0/2  +  1) 

C  =  cos20  tan2  0/2  -  sin20 
T1  =  Z/tan  (0  -  0/2) 
f  ( x )  =  \jx2 C  +  ZBx  -  ZaA 


results  in  the  following  expression: 


f (D/2) 


The  same  equation  yields  V^/4  with  the  substitutions 
B  -B 

D  -D 

T1  -Z/tar(0  +  0/2) 

and  the  total  volume  of  intersection  is  obviously  V  =  4(V  /4  +  V ,lk). 

&  D 

RESULTS 

The  intersection  volume  was  calculated  for  a  fixed  transducer 
spacing  S  and  radius  Al,  for  varying  beam  divergence  angle  ip,  and  beam 
intersection  angle  9  and  the  results  are  plotted  in  Fig.  3.  Inspection 
of  the  plot  indicates  that  the  intersection  volume  differs  from  an  ex¬ 
ponential  function  as  the  beam  intersection  angle  approaches  the  beam 
divergence  angle.  This  result  is  particularly  noticeable  for  the  para¬ 
meter  0  =  6.5°.  The  plots  for  the  two  other  beam  intersection  angles 
are  non- exponential  to  a  lesser  extent  over  the  range  of  values  taken 
for  Ip. 

In  order  to  check  the  numerical  integration  method,  a  Monte  Carlo 
program  was  written  to  calculate  the  intersection  volume  for  0=8° 
and  a  range  of  ip  from  2°  to  4.5°.  For  this  calculation,  the  transducer 


8 


separation  was  fixed  at  3.81  cm  and  the  transducer  radii  at  0.65  cm. 

The  nearly  exponential  dependence  of  the  calculated  volume  as  a 
function  of  if)  is  plotted  in  Fig.  4.  The  plot  shows  the  comparison 
of  the  numerical  integration  and  Monte  Carlo  methods  to  be  nearly 
equivalent  given  the  same  beam  geometry.  The  calculated  intersection 
volumes  found  by  these  two  methods  indicate  that  the  results  of  an 
earlier  study2  are  too  liberal  in  approximating  the  size  of  the  assumed 
elliptical  cross  section. 
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Fig.  4  -  Comparison  of  the  calculated  intersection  volumes  as 
a  function  of  full  beam  width. 
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APPENDIX  I 


This  appendix  Includes  a  copy  of  the  FORTRAN  IV 
coded  program  which  was  used  to  perform  the  intersection 
volume  calculations.  Instructions  for  operating  the  program 
are  included  ir  comment  statements  listed  on  the  printed  output. 
In  addition,  definitions  of  the  input  constants  of  transducer 
size,  intersection  angle,  and  divergence  angle  of  the  acoustic 
beam  are  listed. 
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73/04/13.  il.0C.2G. 

FOG=AM  VOLUME 

00100  prorr-A  '  volv-f(  input, output) 

OOilO*  THIS  CALCULATES  TnK  VOLUME  OF  lNTEr-SECTION  OF 

UOIPO*  THE  PEA*  PATTERN'S  OF  TVP  I DFMTICAL  TF ANSCUCShS.  TnF  LEAN 

(  01 30+  PATTERS  AP.F  Af.pUMED  TO  rZ  COMICAL  ’.  1TH  AN  INCLUDED  AMPLE 

001  /<)*  ''F  rr.r:  EES.  AMI.  THEIR  AXES  IMTEnSIXT  AT  AM  ANGLE  Or 

0(  !*C  +  T-iSTA  DEGREES.  THE  HAD1»»S  OF  T.nE  TRANSDUCERS  Is  At  A\in 

OOlfO*  TIE  5EP  AT.  AT  10  M  OFT  ITEM  TE«MSr*nCV:F  CENTERS  IS 

C  0170+  T.tE  VOLUME  IS  FOUND* PY  C ALCmL ATI  N G  T h K.  AREA  OF  THE 

•  01  +0+  HOM-TONTAL  CROSS-  SECTION  AMP  ''VMEHICALLY  INTEGRA!  I  VP  »nE 

00190*  AREA  *  ITn  S.'FbrECi  TO  THE  2 -COOP  f>  I  MATE  USING  A  nOllERG 

c.  0200+  I  Ml  EFRAVIOM  STir-rtO"TIME  TO  GET  THE  VOLUME. 

00210*  -Op  DERIVATION  OF  THE  FOr.mLAS  AMD  A  MORE  DETAILED 

GCP20*  LESCAI  j-TIOM  OF  THE  CE0ME1RY  SEEi 

v 0230+ 

002*0*  IMr"TSl  rSl  1M  DEGREES 

00250*  THF.TA  IN  DE GHEES 

0 0250*  A 1  IM  AMY  OMITS 

00270*  S  1M  THE  SAME  UNITS  AS  A1 

OOPKQ*  01'TP"TSt  VOL,TME  IM  THE  SAME  .  >1 1 T S  AS  A1 

0  029 C* 

coaoc* 

C  0020  COMMOM  TEMP, C, b» A. D, T2 
C  0330  EMTEEMAL  FCM1 

003 HO  30  PRINT,  *INCLUFED  AMGLF.  OF  TRANSDUCER  BEAM  PATTERN  <PSI>* 

G035G  EFAr. PS11 

0  03  AO  PHIMT*  *AMFLF.  OF  IN?  TEH  SECTION*  OF  BEAMS  (THETA)* 

0  0370  HEAD, THETA 1 

00380  ?FINT,*HADiUS  OF  TRANSDUCERS  (Al)* 

0  0390  HEAD,  A 1 

0  04CG  PRINT,  -DISTANCE  PETVSEM  TFtAMSDI 'CEP.  CEMTEHS  (S)* 

0  0H10  READ,  5 

00420*  CALCULATE  PHI  (THE  AMGDE  Or  TILT  OF  THE  COMES)  AMD  CONVERT 

OCA 30*  ALL  ANGLES  TO  RADIUMS.  FOR  CALCULAT 10M  PURPOSES  PSI/2 

O OH AO*  IS  USED. 

00*50  THETA* THETA 1*3. 14159/iau.O 
00450  PH11*( 1S0.0-THETA1 >/2.0 
0C/70  PH  I *  PHI  1*3. 1  A  I  59/ 180.0 
00480  PSIoPSI 1*3. 1A159/180.0 
0  0490  i  Sl»PSl/2.0 

00500*  CALC-LATE  r  ( THE  DISTANCE  LET?  EEM  VERTICES  OF  TiiS  COMES) 

vOSlu*  CALCULATE  (DISTAMCE  FROM  PLAMF  OF  VERTICES  TO  IMTERSECT  »0 

00520*  OF  AXES. 

00530  X-A1*SIM(THFTA)/TAM(PSI ) 

<•05/0  D-2.0*X+S 
GG55L  £*D*TA.N(  PHI  )/fc.  U 

005A0*  IF  PS1  ♦PH I >90  l.F.GhEF.S  VOLUME  IS  IMFIMITF 

00570  IF(PSI+PHI. CT*  3.141 59/2*  0) GO  TO  1 

00580+  CALCULATE  LIMITS  OF  1'*  .E  LATION  Ult'.P.tZZ)  AML'  COMSTAMTS 

.0  0590*  USED  IM  CALCULATION  (A.F  .C.  TEMr,  Tfe) 

00000  Sl*SIM(rHI) 

COf 10  C  1  =  CQS( PHI ) 

00/20  T=TAM(PSI) 

0  0030  T2«TAM(PHI-PSI ) 

0C./40  A*S1*T*S1*T-C 1*C 1 
00+50  F»2.0*S1*C1*(T*T*1.0) 

OOfAO  C=C 1*T+C1*T-S1*S1 
0  nf-7C  TFv'P-A.0*T*T 
OGGdC  Zi  =  r*T2/2.0 


.  14 

Cep/  (tv-MIt’’  L.uO  'Jotft  tM 

permit  lui’y  rf  production 


0  0690  72=D*S1NC  PHI+FSI  )*SIN'C  rril-Pbl  )/Sl  NC  2.  0*PHI  ) 

0  0700  7  3=  D*TAN  C  PKI+PSI  )/P.O 
0  07  10*  CALCULATE  VA  ( LOVER  VOLUME) 

00720  CALL  XCR0MPC7  1  ,72,*10K-8,  VA.  IERR,  FCN1  ) 

00770*  CHANGE  VARIABLES  FOP.  CALCULATING  VB 

0  07/iC  D=-D 
0  07  60  P---B 

0  07  60  T2=-TA'JC  ,'Hl  +  PSI  ) 

00770+  CALCULATE  OB  CUPPER  VOLUME ) 

00700  CALL  XL  -;OMr  (72.73,1  OE-8,  OR,  l Enft>  FCN  l  ) 

00790*  CALC' .’LATE.  TOTAL  VOLUME  AND  PRINT  RESULTS 

008 00  V=A. 0*( VA+VB) 

C  08  1  0  D=-D 

00820  Pr:I  NT  500,  PS  1  1 ,  TrtET  A  1 ,  rH I  1 ,  S,  A 1 ,  D,  V 

00830  500  FORM  AT  ( *i'SI  =  *,F1G.2.+  TRETA=  *,F10.2,+  rrii  =  +  .F10.2, 

0  08  A0  +  /  *  S=  *  ,  F  1 0.  2,  *  A 1 =  +.F1Q.2,*  D=  +.F1G.2/ 

00860+  * VOLUME  OF  INTERSECTION  =  +.E15.8) 

0  08  60  GO  TO  20 

00870  1  PT I NT  501, PS1 1,  THETA1, PHI  1,  S,  Al,  D 

00880  501  FO.-.MATC  *  PSI  =  *,F10.2,+  THETA=  *,F10.2,=  PHI=  *,F10.2/ 
00890+  *S=  *,F10.2,+  A 1 =  *,F10.2,*  D=  *,F10.2/ 

009U0+  * VOL "ME  OF  INTERSECTION  IS  INFINITE*) 

00910  20  PRINT  50J 
0  0920  505  FORM  AT (  //  ) 

0 09 30  GO  TO  30 
009 AO  END 

00950  FUNCTION  FCN  1(7  ) 

00960*  THIS  FUNCTION  CALCULATES  THE  AREA  OF  CROSS-SECTION  AT 

00970*  HEIGHT  7.  IT  FIRST  CALCI'LATES  THF  LIMITS  AT  l, RICH  THE 

0  0980*  AREA  INTEGRAL  MUST  BE  EVALUATED  AND  THEN  CALLS  FCN'P  AT 

00990*  THESE.  LIMITS  TO  MAKE  THE  EVALUATIONS. 

OlOCO  COMMON  TEMP, C, t, A, D, T2 
01010  Y 1 =E/2. 0 
0  1020  Y  2=7  / T2 

0  1030  Y1=C*Y1*Y1+7*Y1*P+7*7*A 
0  1  O/iO  Y2  =  C*Y2*Y2+7*Y2*R  +  7*7*A 
01050  I F ( AfcS ( Y 1 ). LI . ICE- 10) Y 1=0.0 
0  1060  IF (  ABSC Y2).LT. 1 OE- 1 0 ) Y2=U. 0 
01070  IF<Y1.LT.0.0)PRIN1,*Y1»2=+, Y1 
0  108  0  IFCY2.LT .0.0) PRINT, *Y2» 2=*, Y2 
0  1 090  Y 1  =  SORT ( ABSC  Y 1 ) ) 

0  1 100  Y 2=  SORT C  AF  S  C  Y2 ) ) 

OHIO  FCN  1  =  FCN2  C  Y2,  7  )  -  FCN2C  Y  1 , 7  ) 

01120  RETURN 
01130  END 

0  1  1  AC  FUNCTION  FCtNR.C  X,  7  ) 

01150*  THIS  FUNCTION  EVALUATES  THF  AREA  INTEGRAL  AT  X,7. 

01160  COMMON  TEMP, C,F. A, D, T2 
0  1170  HA P=7  *  7  *  TEMP/ C -* . 0  +  C ) 

0  1180  F=X/SOP.T  C  RAD) 

0  1190  F2=F.AD-X*X 

0  1200  1FCF.LE.  1.0)  GO  TO  1 

0  1  2  1  u  PRINT,  *AHC-  OF  AS  IN  >  1*,F 

01220  F= 1 • 0 

01230  v2=0« 0 

0  1  2 AO  1  “’CM 2=  (  ,\*SL:.T  C  F2  ) -RAT +ASINC  F)  )/<  SORT  (  -0*2) 

0  1250  RETURN 
C  1260  END 
h  EAl  Y . 
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